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Abstract 

Qfs ' In a previous work regarding the interaction of two p{770) resonances, the /2(1270) ( J^*^ — 2++) resonance 

was obtained dynamically as a two-p molecule with a very strong binding energy, 135 MeV per p particle. 
In the present work we use the pp interaction in spin 2 and isospin channel to show that the resonances 
Q_^. P3(1690) (3~), /4(2050) (4++), P5(2350) (5"-) and /6(2510) (6++) are basically molecules of increasing 

number of p(770) particles. We use the fixed center approximation of the Faddeev equations to write the 
multi-body interaction in terms of the two-body scattering amplitudes. We find the masses of the states 
(-H ' very close to the experimental values and we get an increasing value of the binding energy per p as the 

number of p mesons is increased. 
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INTRODUCTION 



In the last decade, the chiral unitary approach has shown that many hadronic resonances 
can be obtained dynamicahy from the interaction of hadrons. This has been done through the 
implementation of unitarity in coupled channels using a lowest order chiral Lagrangian, to the 
point that these resonances can be interpreted as meson-meson or meson-baryon molecules [ll4lO|. 



and it has shed new light into the issue of the nature of the scalar mesons, among others. The 
interaction of pseudoscalar mesons among themselves and meson-baryon interaction has given way 



recently to the interaction of vector mesons among themselves Ill4l3l| and the interaction of vector 






mesons with baryons [IJ, ll5| , where the interaction is evaluated within the techniques of the chiral 



unitary approach starting from a lowest order hidden gauge symmetry Lagrangian 16l4l9l|. In the 



work [ll[ it was found that the interaction of two p{770) mesons in isospin 7 = and spin 5 = 2 
was strong enough to bind the pp system into the /2(1270) (J^*-" = 2++) resonance. The nature of 
this resonance as a />(770)p(770) molecule has passed the tests of radiative decay into 77 120], the 



decay of J/^ into (^{(p) and /2(1270) (together with other resonances generated in [I3]) [211], and 



J/\I' into 7 and /2(1270) (and the other resonances of |l2l]) [22 

The /2(1270) obtained in ref. [11] as a pp quasibound state or molecule implies a very large 
binding energy per p meson, about 135 MeV. This occurs only for spin S=2, where the two spins 
of the p are aligned in the same direction. In view of this strong pp interaction, some natural 
questions arise: i) is it possible to obtain bound systems with increasing number of p mesons as 
building blocks? These systems with many /o's with their spins aligned in the same direction would 
make a condensate, with features similar to a ferromagnet; ii) If so, is there a limit in the number 
of /o's or, even more interesting, the mass of the multi-/? system saturates at some number of p 
mesons?. In this latter case then it would be energetically "free" to introduce new p mesons into 
the system. 

The condensates made out of mesons have been advocated some times, and concretely, the issue 
of pion condensates was popular for some time [23!] and kaon condensation has also attracted much 
attention J24]. 

Regarding the question i), in the PDG [25( there are intriguing mesons with large spin, of the 
p and /o type, whose quantum numbers match systems made with 3, 4, 5 and 6 p mesons with 
their spins aligned. These are the />3(1690) (3~), /4(2050) (4++), /95(2350) (5~) and /6(2510) 
(G^"*") resonances. If these resonances were essentially multi-/? meson molecules, they would have 
a binding energy per p of about 210, 260, 305 and 355 MeV, respectively. This increasing value as 
more p's are added to the system connects with question ii). 

The main aim of the present work is to address these questions. Technically, we would have to 
solve the Faddeev equations [2Q] for a state of three p's to start with. Systems of three mesons. 



concretely KK(j), have been addressed recently with Faddeev equations 27|]. States of two mesons 



and one baryon have also received recent attention 2814311]. In the three p mesons state, the 



fact that the two p meson system with S=2 is so bound, makes it advisable to use the fixed center 



approximation to the Faddeev equations (FCA) 32H35I] in order to obtain the scattering amplitudes 



of one p with the /2(1270) state. The FCA requires the knowledge of the wave function of the 
bound state of the target. This information can be obtained using a recent method that connects, 
in an easy and practical way, the wave functions with the scattering matrices of the chiral unitary 
approach [36] . Proceeding by iterations we build up states with an extra p meson starting from the 
former state. In this way the multi-p resonances are generated, which show up as prominent bumps 
in the different scattering amplitudes. The iterative method suggest a way to extrapolate to many 
p states and we develop an analytical method to evaluate the mass of the n-p state for very large 
n if only the single-scattering contribution is considered in the Faddeev equations, which should 
be viewed only as suggestive of what might happen in the limit of a large number of p mesons. 



II. UNITARIZED pp INTERACTION 

The most important ingredient in the calculation of the multi-p scattering is the two-p interac- 



tion. In this section we briefly summarize the model of ref. 11[ to obtain the unitary pp scattering 
amplitude. (We refer to [11] for further details). 



The pp potential is obtained from the hidden gauge symmetry Lagrangian 16l4l9l| for vector 
mesons, which, up to three and four vector fields, reads: 

£(3^) = igiiV^'d.V^ - d.V^vnvn , (2) 

where Vu is the SU{3) matrix containing the vector-meson fields and the coupling constant g is 
g = Mv/2f with / = 93 MeV, the pion decay constant. The Lagrangian of Eq. ([1]) gives rise to 
a four vector meson contact term and that of Eq. (I5|) to a four vector meson interaction through 
the exchange of an intermediate vector meson in the t and u channels (the s-channel gives rise to 
a p-wave that we do not consider, only the important s-wave part is studied). 

From these Lagrangians, a potential V can be obtained to which the contact and p-exchange 
terms contribute. For the present work only the spin S = 2 and isospin 1 = 0,1 = 2, are necessary: 

y (7=0,5=2) (^) = -V - 8g' (^ - l) ~ -20ff2 

where the last terms are the approximate values at threshold in order to give an idea of the weight 
and sign of the interaction. The pp S = 2, I = is strongly attractive. This is the most important 
reason to obtain a bound pp state with these quantum numbers as we explain below. 



Further contributions to the previous potential where considered in ref. Ul], out of which only 



the box diagram, which accounts for the two-pion decay mode, was relevant, and only for the 



imaginary part of the potential. Explicit expressions can be found in ref. [ll[ |. 

With this potential the total pp scattering amplitude can be obtained. In order to extend 
the range of applicability of the interaction to the resonance region, the implementation of exact 
unitarity is mandatory. In this case, we use the Bethe-Salpeter equation where the kernel is the 
potential V described above: 

- '^ (4) 



for each spin-isospin channel. In Eq. HI G is the unitary bubble or the pp loop function [2, y] 

^^""^ " V {2TT)'^p^-mj + ie{Q-pf-mj + ie ' ^^^ 

with Q = (^/s,0). This loop function can be regularized by means of dimensional regularization 
or using a three-momentum cutoff, Pmax = A: 



of + mjy/^. 

In order to consider the width of the p particles inside the loop, a convolution with the two p 
meson spectral functions is done to Eq. (j6l): 



where Wj = (p ? + m^ 



^ Mmp+2ToN? Amp+2ToN? 

Gis) = Yn '^^1/ ds2G{s,^/s^,^/s^)Sp{sl)Sp{s2) (7) 

■'V J(m„-2rnArV J(m„~2rnAr)^ 



^ Mmp+2ToN? Amp+2ToN) 

/ dsi / 

'(mp-2roiv)2 J{mp~2ToN? 

where Sp{si) is the p meson spectral function 

Sp{si) = Im 

IT 



Si-mj + impTp{y/s'i) 
with Tp{^/s) the p-meson energy dependent width 



(8) 



/ s — Am \ 

r,{V-s) = roN -^-±2) (9) 



^?7ip — Am^ J 
and Ton the on-shell p meson width. In Eq. ([7|) A/" is a normalization factor given by 

f{mp+2VoN? 

M= / dsSp{s) (10) 

J{mp-2roiv)2 

The cutoff A is the only free parameter in the whole model and is chosen such as to produce 
the peak of |rp at the experimental mass of the /2(1270). This implies A ~ 875 MeV, which is of 
a natural size [81], about 1 GeV. 

In fig. [H the modulus squared of the 5 = 2, / = scattering amplitude, T^^^^'^^"^', is plotted. 
The resonance structure of the /2(1270) resonance is clearly visible. 

III. MULTI-BODY INTERACTION 

We are going to use the fixed center approximation of the Faddeev equations in order to obtain 
the interaction of a number of p mesons larger than two. 

We will illustrate the process for the interaction of three mesons and will give the expression 
obtained analogously for other number of mesons. For the three p system, we will consider that 
two of the p mesons are clusterized forming an /2(1270) resonance, given the strong binding of the 
72(1270) system. This allows us to use the FCA to the Faddeev equations. 

The FCA to Faddeev equations is depicted diagrammatically in fig. [2j The external particle, 
the p in this case, interacts successively with the other two p mesons which form the pp cluster. 
The FCA equations are written in terms of two partition functions Ti, T2, which sum up to the 
total scattering matrix, T, and read 

Ti = ti + tiGoT2 

T2 = t2 + t2GQTi 

T = Ti + T2 (11) 
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FIG. 1: Modulus squared of the pp scattering amplitude with total spin S = 2 and isospin / = 



where T is the total scattering amplitude we are looking for, Tj accounts for all the diagrams 
starting with the interaction of the external particle with particle i of the compound system and 
ti represent the pp unitarized scattering amplitude of a p'^ with any of the other p in the I = pp 
system. The schematic representation is depicted in fig. [2j 

Fig. [2^) represents the single-scattering contribution and fig. [2b) the double-scattering. The 
contributions of fig. [2^ and b are the two first contributions of the Faddeev equations. 

In the present case, since both 1 and 2 are p mesons we have Ti = T2 an thus the system of 
equations is just reduced to a single equation 



T 



h + tiGoTi 
2Ti 



(12) 



A. Single-scattering contribution 



The amplitude corresponding to the single-scattering contribution of fig. [2^ comes just from the 
ti term of Eq. ^,T = 2ti. 

In order to write this expression in terms of the 1 = and 1 = 2 unitarized amplitudes {tpp~ , 
tpp ) of Eq. ([!]), let us consider a cluster of two p mesons in isospin / = 0, the constituents of 
which we call mesons 1 and 2. The other p meson will be meson number 3. The two p mesons 
forming the /2 are in an / = state 



\pp)i=o = --^IpV + p-/^+ + pV) = 7f ( l(i'-i)) + K-i'i)) - l(O'O)) 



(13) 



where the kets in the last member indicate the Iz components of the 1 and 2 particles, \{Iz ,Iz ))• 

(3) 
We take the p meson number 3 in the state |(/^ )) 

|p+) = -|(+l))- (14) 
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FIG. 2: Diagrammatic representation of the fixed center approximation to the Faddeev equations. Diagrams 
a) and b) represent the single and double scattering contributions respectively. 



The scattering potential in terms of the two body potentials V31, V32 is: 



-((+1)1 i= (((+1, -1) + (-1,+1) - (0,0)1) ) (^31 + ^32) 
■|(+1)) -^ (|(+1, -1) + (-!,+!) - (0,0))) 
^(((2 +2), -1) + (-^(20) + i=(10) + i=(00), 1) - (i=(2 +1) + i=(l +1), 0) 



V, 



31 



((2 +2), -1) + (-L(20) + i=(10) + i=(00), 1) 



:^(2+l) + ^(l+l),0) 



+ WtE^''^ + 7=2^''^ + 71^'"^' '^ + ^^' +'^' -'^ - ^T2^' "-'^ + 7=2^' ^'^' "^ 



V32 



((-^(20) + -^(10) + i=(00), 1) + (2 +2), -1) - (-L(2 +1) + -L(l +l),0)y (15) 

where the notation followed in the last term for the states is ((1*°*^^/*°*^', /^)|T^j|), where I^"^'^^ 



means the total isospin of the ij system and k ^ i,j (the spectator p). 

This leads, in terms of the 1 = and 1 = 2 unitarized amplitudes {tpp~ , tpp~ ), to the following 
amplitude for the single scattering contribution: 



^(5tr^+t(rT (16) 



2 

ti = - 
9 

where we have added an extra 2 factor in order to match the unitary normalization of ref. [Uj in 



pp- 



It is worth noting that the argument of the function T{s) is the total invariant mass energy s, 
while the argument of ti is s', where s' is the invariant mass of the p meson with momentum k 
and the p meson inside the /2 resonance with momentum pi and is given by 

s' = {k + p,f = ]^{s + ^ml-Ml) (17) 

For latter applications, let us write the general expression of s' for the interaction of a particle 
A with a molecule B with n equal building blocks h. Then, s' represents the invariant mass of the 
particle A and a particle h of the B molecule and is given by 

s' = -(s-Ml- Ml) +Ml + ml (18) 

n 

where M^/^) is the mass of the A{B) system and ttt-^ is the mass of every building block of the B 
molecule. 

Let us consider the wavefunctions of the incident and outgoing p particles being plane waves 
normalized inside a box of volume V and let us call tpi the wavefunctions of the p mesons inside 
the /2 resonance. The S'-matrix for the process of Fig. [2^ is written as 

which we can multiply by the identity 

cPx'^2{x')Mx') = 1- (20) 

The integration of the time component x^ provides the energy conservation at the interaction 
point x: 

I dx%-'P°^\'P''''\-'''°'=\"''°'=° = 27r Sip'i + A:0 - p'? - k"^) = 2tt ^(/fcO + Ef, - k"^ - E'j^), (21) 

tiere in the last step we have assumed P2 = P2 ) 
spectator (impulse approximation). We can take 



where in the last step we have assumed P2 = P2 ) ^^ corresponds to having the second particle as 



Mx)Mx') = -^e'^f^-^^f,{f), (22) 

vV 



with ^j2 the wave function of the /2(1270) cluster and 



f = x — X . (23) 



and then we get for the spatial integrals 

f d'R e^^/2-^e-''^/2- V-^e-*^^'-« = {2^f5{k + Kf, - k' - ^J (24) 

and 

jd\^f,{r)^f,{r)e^He~^^^'-i = Ff,['^) ^ FjM = 1 (25) 

where Fj^ is the /2(1270) form factor normalized to unity neglecting the k, k' momenta, which we 
take equal. 

Hence the S'-matrix for the single scattering term is given by 

S^'^ = -^*ii^^^i=7i=^(2vr)^ 5{k + Kj, - k'' - K'f^). (26) 

and recall we must sum ti + ^2 — ^ 2ti. 

B. Double-scattering and resummation contribution 

We are going to evaluate the amplitude of the double-scattering contribution (fig. [2)3) in a 



similar way as in the case of the kaon deuteron interaction in 35|, |37 1 . 
The S-matrix can be written as 

5(^) = [d'x /d^x'^e--?-Vi(x)-^e-'?-Vi(xl^e--^-'V2(x0^e-'^-'V2(x') 
J J y^2u}p^ V^^p'i V2a;p2 V^^p'2 

^^^-L^^^, / ^ \ / , (-.t,)(-iti) (27) 



^J%J^ ^/2^y J {2tt)^ q^ - mj + ie 

The integrations of the time components x^ and x' provide the energy conservation at the two 
interaction points x and x': 

I dx^e-'P°^°e'P'°''°e-'^°^°e'^°^° = 27r 5{p^, + k^ - p'\ - q') 

I dx'°e-^°-'"e^P'°-'°e^^'°-'°e-^"-'° = 2vr 6{pl + q' - A " A:'°) (28) 

We implement now the change of variables (x, x') — )• {R,r) of Eq. (I23p . The R integral gives 
the same expression as in Eq. (j24p , and the r integral gives rise to 

d\e-^(^^~^y^^^f2iff ^ Fj,[q- ^) (29) 

where Fj^ [q — {k + k')/2] is the /2(1270) form factor introduced above. 



The final expression for tlie 5-inatrix for the double scattering process is 

1111 1 1 1 



^ ^' V^ V^ y^ V^i V^i V^2 V^P 



and we will take q^ at the /2 rest frame, q^ = {s — rm? — M? )/{2Mf^), where we have considered 
pi^ = p'^ and p2^ = p'2 which is true in average. In Eq. (|30p we have also taken into account that 
{k + k')/2 = in average. 

For the evaluation of the form factor of the /2 resonance we follow the approach of [So] . In this 
work it is shown that the use of a separable potential in momentum space of the type 

V = ve{A - q)e{A - q') (31) 

where A is the cutoff used in the theory for the scattering of two particles and q, q' are the modulus 
of the momenta, leads to the same on shell prescription for the scattering matrix as is used in the 
chiral unitary approach. The on shell prescription converts the coupled integral equations for the 
scattering matrix into algebraic equations, and similarly, the wave functions can be easily obtained 
in terms of an integral. The wave function in momentum space is written as 

(P IV) = V ®^^-P^ I d\{m, (32) 

E-ujp{pi)-u:p{p2) Jk<A 



where u!p{p) = \/p'^ + rn'l, and in coordinate space as 



i^\i^)=lj^^''-'ip\^)- (33) 

The final expression for the form factor of Eq. (|29|) is then given by 



where the normalization factor M is 



1 



M= / d-^p ^. (35) 

Jp<A {Mf,-2up(p)f 

In fig. [3] we show the form factor of the /2 resonance. The condition \p — (f\ < A implies that 
the form factor is exactly zero for q > 2A. Therefore the d\(f\ integration in Eq. pop has an upper 
limit of 2 A. 

We must now face the issue of normalization in our formalism. We use Mandl-Shaw [38] nor- 
malization for the fields and hence the S'-matrix for pf2 scattering is written as 



' = -•^'''<'>v^^ vSTTSI^^'^''*'** + ">' - *'° - "'"^ <^«' 
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FIG. 3: Form factor of the /2(1270) resonance 
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FIG. 4: Real and imaginary parts of the Go function, Eq. ([38|) 



but this should be compared with expressions Eq. ()26|) for the single scattering and Eq. (|30p for 
double scattering. Summing the two partitions T\ and T^ we find that 



^ph 



4(ti+titiGo) 



(37) 



where we have made the assumption that in the /2 rest frame, where we evaluate the amplitude, 
Ibjp ~ Mj2 , and Go is given by 



Go 



1 



^^h 



(2^ 



FfM)- 



1 



,02 



^2 — 771,2 _|_ ^g 



(38) 



In fig. m we show the real and imaginary parts of the Go function. Note that close to the 
threshold it has the typical shape of a two meson loop function, pf^ in this case, but it is smoothed 
towards zero at higher energies due to the form factor. 



10 



Equation (j37p represents the two first terms of the series expansion of 4fi/(l — tiGo). Actually, 
if we consider further number of scatterings in the expansion of Tpf^ of the FCA, (see diagrams d 
in fig. [2]), we get 



^p/2 — 4(^1 + ^iG'o^i + tiGotiG^ti + tiGotiGotiGoti + 



4ti 



tl\s') 



d?q 
where s' is given in Eq. ()17p . 



FfM- 



,0 _ g*2 _ ^ _|_ ^g 



1 - Goii t^i - Go 



(39) 



C. Larger number of p mesons 

For the interaction of up to 6 p mesons we can follow a similar procedure as in the previous 
subsections but considering the interaction of two different clusters. For the interaction of 4 p 
meson we can calculate the interaction of 2/2(1270) resonances given the strong tendency of two p 
mesons to clusterize into an /2. Advancing some results that we will show later, this four p state 
gives rise to the /4 resonance. Thus, analogously, for 5/7 we can consider the interaction of one p 
meson with an /4. And for Qp we can consider the interaction of an /4 with an /2. 

Therefore the amplitude for the interaction of a cluster A with a cluster B made of two equal 
components b is given by 



where 



t{s- A,B)=A [t-\s'{s- A, b); A, b) - Go(s; A, B)] ' 



Go{s-A,B) = ^ f^F{q;B) 



Mb J (27r)2 



q'^{s;A,Bf -q^-Ml + ie' 



(40) 



(41) 



F{q,B) 



., d'p- 



^ ■'\/-f<A' Mb - 2 Jp2 +mlMB- 2J\p- e + ml 



(42) 






(43) 



and, from Eq. (jlSp. 



q"is;A,B) 



Ml 



Ml 



2Mb 



(44) 



s'{s- A,b) = -{s-Ml- Ml) +Ml + ml 



(45) 



Note that it is not necessary that the cutoff A' be the same in all the cases as the A used for 
the /2 case. The cutoff A used in Eq. Q for the pp loop function, which is the same appearing in 
the momentum integral to get the /2 form factor in Eq. ^^ [3^], can be interpreted as the typical 



11 



maximum momentum that each p can reach inside the /2 molecule. For the /4 case we can argue 
that the maximum value would be like in the /2 case but scaled by the typical momentum of the 
/2 components inside the fi molecule. The typical three- momentum of the components of the /4, 
74, is of the order of 



14 ^ \l — + Mf^B ; B = Mf^- 2M/2 (46) 



where B is the binding energy of the f^ and analogously for the f2'- 



4 
This gives for the cutoff of the /4 



A' 



^2 

l2^\l-r + mpB ■ B = Mf^-2mp (47) 



A^ 



~ 1500 MeV. (48) 

u " 

While this is just a very rough estimation, this gives us an idea of the order of A'. In any case this 
only affects the evaluation of the 5p and 6/9 system. In the numerical evaluation we will consider the 
range A' ~ 875 — 1500 MeV to have an idea of the uncertainties from this source. But, advancing 
some results, the dependence of the mass of the systems with this cutoff is small. 

D. Arbitrary number of p mesons in single scattering approximation 

For and arbitrary number of p mesons, it is possible to obtain a simple analytic expression for 
the mass of the multi-p system if only the single scattering mechanism is considered, which is the 
first order mechanism. Of course this is just a toy approximation since, as we will see in the results 
section, the multiple scattering is important, but it serves to make some interesting qualitative 
arguments. 

In the single scattering approximation, Eq. (j40p takes the form 

t{s;A,B)^At{s'{s;A,b);A,b) (49) 

That means that, for instance, the amplitude tpj^ is just proportional to the tpp amplitude but 
evaluated at an energy s' shifted with respect to s due to the fact that one of the p's involved in the 
pp scattering is bound into an /2 system. In general, the interaction amplitude of a number Up of 
p mesons is proportional to the pp amplitude with an energy obtained considering that one of the 
p mesons is bound into an {up — 1) molecule. Therefore one can obtain recursively the amplitude 
for the Hp system. Because of that, the shape of |t(,s; A, B)'\^ is the same as that of |ipp(s)p but at 
a shifted energy. The s value at which |ipp(s)p has the maximum is precisely M? . The value of 
s appearing in i(s; A, B) of Eq. ([l9j) is the value that we can assign to the mass of the Up system, 
M{np). Therefore, applying recursively the above condition one can obtain a general expression 
for M{np) in the single scattering approximation: 

M{npf = ^np {rip - 1) M\ - Up (up - 2) m^. (50) 

We can also define a binding energy per p as 

E(np) = "P"^P - ^("P) (51) 

Tin 



which will be used for later discussions. 
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FIG. 5: Modulus squared of the unitarized multi-p amplitudes. Solid line: full model A'| , — 1500 MeV; 
dashed line: full model A'L = 875 MeV; dotted line: only single-scattering contribution. (The dashed and 
dotted lines have been normalized to the peak of the solid line for the sake of comparison of the position of 
the maxima) 



IV. RESULTS 



In fig. [5] we show the modulus squared of the amphtudes for different number of p mesons 
considering only the single scattering mechanisms (dotted line) and the full model (solid and dashed 
lines) . The difference between the solid and dashed lines is the value of A' L of Eq. (j38|) needed in 
the evaluation of the 5/9 and 6/9 meson systems (1500 MeV in the solid line, 875 MeV in the dashed 
one). The dotted and dashed curves have been normalized to the peaks of the corresponding 
full result for the sake of comparison of the position of the maximum. The difference between 
the dashed and solid lines can be considered as an estimate of the error but the variation in the 
position of the maximum is small. 

We clearly see that the amplitudes show pronounced bumps which we associate to the resonances 
labeled in the figures. The position of the the maxima can be associated to the masses of the 
corresponding resonances. 

In table H] the values of the masses of our generated multi-/> systems are shown in comparison 
with the experimental values at the PDG [25]. The two values for the p^ and /e masses in the full 
model column correspond to the different values in A'L. as explained above. In the last column 



the binding energy per p meson, E{np) 



[n„mn 



M{np))/np is also shown. 



In fig. [6] we show graphically the results for the masses of table HI 
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np 




mass, PDG [25] 


mass, only single scatt. 


m.ass, full m.odel 


E{np) 


2 


/2(1270) 


1275 ±1 


1275 


1285 


133 


3 


P3(1690) 


1689 ± 2 


1753 


1698 


209 


4 


/4(2050) 


2018 ±11 


2224 


2051 


263 


5 


P5(2350) 


2330 ± 35 


2690 


2330-2366 


302-309 


6 


/6(2510) 


2465 ± 50 


3155 


2607-2633 


337-341 



TABLE I: Results for the masses of the dynamically generated states. 
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FIG. 6: Masses of the dynamically generated states as a function of the number of constituent p(770) mesons. 
Up. Only single scattering contribution (dotted line); full model (solid line); experimental values from the 
PDG[25], (circles). 



We can see from the results that the single scattering mechanism produce quahtatively the 
resonances but the positions of the masses do not agree with the experimental values by differences 
ranging from about 60 MeV for the /O3 to 700 MeV for the /g. The situation is drastically improved 
when the multiple scattering is considered. In this case, the agreement with the experimental 
values of the masses is remarkable. Quantitatively, the full model is essentially compatible with 
the experimental values within errors except for the f% where the discrepancy is about 150 MeV, 
which is still quite remarkable, given the high mass and width of the resonance. The typical 
discrepancy with the experimental masses is of the order of 1%, (5% for the /e). 

It is worth stressing the simplicity of our approach and the absence of parameters fitted in the 
model. To be more precise, only the value A = 875 MeV of the cutoff of the pp loop function was 
chosen in ref. HI] to agree with the experimental f^ pole position. No further adjustments have 
been done in the present work. 

In principle, the widths of the bumps can be associated to the the widths of the resonances if 
they were Breit-Wigner like shapes, which is clearly not the case. This means that the amplitudes 
contain much non-resonant background which our model generates implicitly through the non- 
linear dynamics involved in the unitarization procedure. That means that the extraction of the 
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widths of the resonances from our amplitudes is just very quahtative: 200, 350, 900 and 1500 MeV 
for /93, /4, P5 and /g respectively. The order of magnitude agree with the experimental value of 
the PDG [2S], 161 ± 10, 237 ± 18, 400 ± 100 and 255 ± 40 respectively, except for the two heaviest 
states. However, it is worth noting that, by looking at fig. O in these heaviest states much of the 
strength of the amplitude off the peak could be interpreted as a background, as would be the case 
in an experimental analysis of a distribution like the one obtained in fig. [5l in which case the actual 
width of the resonance would be significantly reduced. 

Let us address again the problem of an arbitrary large number of p mesons. A natural question 
looking at fig. [6] is if the curve of the masses saturates for a large enough number of p mesons. 
That would imply that it would be energetically free to add an extra p meson to the system. For 
the single scattering case we can analyze the problem with the help of Eq. (I50p . The saturation 
would occur for 

ml - M? /4 
sat ml - M2 /2 ^ ^ 

which never happens for the actual p and /2 masses, since Up of Eq. (152p gives a negative value. 

sat 

However it is worth noting in fig. [6] that the single scattering is just a bound limit and that the 
multiple scatterings tend to decrease importantly the mass of an np system. If such decrease is 
enough to eventually reach the saturation condition cannot be answered with certainty within the 
present model since we do not go beyond Up = 6. However, it is worth noting the large value of 
the binding energy per p, see last column of table HI Already for Up = 6, it is almost half the value 
of the p meson mass. That means that the creation of the 6 p meson system gives back half the 
mass of all the particles involved which is quite a lot of energy. 

The binding energy per p evaluated using the single scattering approximation, Eq. (|5ip . tends 
asymptotically to 



M2 
lim E(np) = mp - \ -^^ - ml c^ 315 MeV. (53) 

rip— j-cxD y 2 

However this value is already reached at np = 6 if the multiple scattering mechanisms are consid- 
ered. 

If the pp interaction were a little bit stronger, such that Mf^ ~ \/2m,p = 1096 MeV, then the 
saturation would be reached already considering only the single scattering. And, in order to get 
saturation for Up = 6 with only single scattering, the mass of the /2 resonance should be just 
slightly smaller, ~ y^20/llmp = 1056 MeV. Of course this is just a qualitative reasoning since the 
width of the system would eventually increase with the number of p mesons, making the system 
fade away rapidly. The former discussion is obviously rough and speculative of what might happen 
for large np systems. What remains as quantitative results from the present study is the fact that 
the /2(1270), /93(1690), /4(2050), /95(2350) and /6(2510) can be essentially considered as multi-p 
molecules with increasing number of p mesons. 

V. CONCLUSIONS 

In the present work we claim for the first time that the p3(1690), /4(2050), /95(2350) and /6(2510) 
resonances can be interpreted as multi-p states of 3, 4, 5 and 6 p mesons respectively, with their 
spins aligned. The main idea stems from the fact that in ref. [11| it was found that the interaction 
of two p{770) mesons in isospin 1 = and spin S" = 2 is very strong, to the point to bind the two p 
mesons forming the /2(1270) resonance. This elementary pp interaction is obtained implementing 
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unitarity, using the techniques of the chiral unitary approach, with a potential obtained from a 
hidden gauge symmetry Lagrangian for the interaction of two vector mesons. For the multi-p 
systems we evaluate the scattering amplitudes for the interactions of two clusters made up of p- 
mesons. To this purpose we use the fixed center approximation to the Faddeev equations which 
considers the multiple scattering steps in addition to the single process where each p meson interacts 
with all the rest of p mesons within the cluster. 

The position of the maximum in the modulus squared of the amplitudes can be associated 
with the masses of the corresponding resonances. It is worth noting that the model has no free 



parameters once a cutoff is chosen in ref. Ul to obtain the experimental mass of the /2(1270) 
resonance. 

The values of the masses that we obtain are in very good agreement with the experimental 
values of the masses of the resonances considered in the present work, the p3(1690), /4(2050), 
yO5(2350) and /6(2510). This is a remarkable fact given the simplicity of the underlying idea. 

The states obtained have an increasing binding energy per particle which induces one to spec- 
ulate on the possibility that for a given number of p mesons it would cost no energy to produce 
a new p meson inside the meson condensate state. However, simultaneusly we observe that the 
width of the new systems also increases with the number of p mesons, to the point that for Up = Q 
the width is already very large. It might as well be that one has reached an experimental threshold 
and that new multi p states, that in principle could be created, have such a large width that they 
escape present detection techniques. In any case, the claims made here that the already observed 
states up to J = 6 correspond to multi p states is a novel idea worth consideration. New studies 
with different formalisms and different points of view would be most welcome, as well as possible 
experimental tests which could help unveil the real nature of these states. 
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